#Plot parameters
fig_ndvi <- ggplot(aes(x = Plot_rev, y = NDVI, fill = Plot_rev), data = Plot_variables_s) +
#geom_violin(trim = FALSE, alpha = .25, colour = NA) +
geom_jitter(width = .15, aes(color = Plot_rev), alpha = .5, shape = 16, size = 3) +
stat_summary(fun.data = mean_cl_boot, geom = "errorbar", width = .15, color = "grey25") +
stat_summary(fun.y = mean, geom = "line", size = .5, aes(group = Location), color = "grey25") +
stat_summary(fun = mean, geom = "point", size = 4, aes(fill = Plot_rev), shape = 22) +
scale_color_manual(values=c("#3aba1e","#a1da27","#f0f123","#f5a527","#db2f14")) +
scale_fill_manual(values=c("#3aba1e","#a1da27","#f0f123","#f5a527","#db2f14")) +
ylim(0,1000) +
theme_bw() +
theme(strip.background=element_blank()) +
labs(x ="Grazing pressure", y ="Primary production (NDVI)") +
facet_grid(~Location)
fig_biomass <- ggplot(aes(x = Plot_rev, y = Biomass, fill = Plot_rev), data = Plot_variables_s) +
#geom_violin(trim = FALSE, alpha = .25, colour = NA) +
geom_jitter(width = .15, aes(color = Plot_rev), alpha = .5, shape = 16, size = 3) +
stat_summary(fun.data = mean_cl_boot, geom = "errorbar", width = .15, color = "grey25") +
stat_summary(fun.y = mean, geom = "line", size = .5, aes(group = Location), color = "grey25") +
stat_summary(fun = mean, geom = "point", size = 4, aes(fill = Plot_rev), shape = 22) +
scale_color_manual(values=c("#3aba1e","#a1da27","#f0f123","#f5a527","#db2f14")) +
scale_fill_manual(values=c("#3aba1e","#a1da27","#f0f123","#f5a527","#db2f14")) +
ylim(0,300) +
theme_bw() +
theme(strip.background=element_blank()) +
labs(x ="Grazing pressure", y ="Primary production (biomass)") +
facet_grid(~Location)
fig_cover <- ggplot(aes(x = Plot_rev, y = Plant_cover, fill = Plot_rev), data = Plot_variables_s) +
#geom_violin(trim = FALSE, alpha = .25, colour = NA) +
geom_jitter(width = .15, aes(color = Plot_rev), alpha = .5, shape = 16, size = 3) +
stat_summary(fun.data = mean_cl_boot, geom = "errorbar", width = .15, color = "grey25") +
stat_summary(fun.y = mean, geom = "line", size = .5, aes(group = Location), color = "grey25") +
stat_summary(fun = mean, geom = "point", size = 4, aes(fill = Plot_rev), shape = 22) +
scale_color_manual(values=c("#3aba1e","#a1da27","#f0f123","#f5a527","#db2f14")) +
scale_fill_manual(values=c("#3aba1e","#a1da27","#f0f123","#f5a527","#db2f14")) +
ylim(0,100) +
theme_bw() +
theme(strip.background=element_blank()) +
labs(x ="Grazing pressure", y ="Primary production (plant cover)") +
facet_grid(~Location)
#Print
grid.arrange(fig_ndvi, fig_biomass, fig_cover, nrow = 3)
fig_rich <- ggplot(aes(x = Plot_rev, y = Richness, fill = Plot_rev), data = Plot_variables_s) +
#geom_violin(trim = FALSE, alpha = .25, colour = NA) +
geom_jitter(width = .15, aes(color = Plot_rev), alpha = .5, shape = 16, size = 3) +
stat_summary(fun.data = mean_cl_boot, geom = "errorbar", width = .15, color = "grey25") +
stat_summary(fun.y = mean, geom = "line", size = .5, aes(group = Location), color = "grey25") +
stat_summary(fun = mean, geom = "point", size = 4, aes(fill = Plot_rev), shape = 22) +
scale_color_manual(values=c("#3aba1e","#a1da27","#f0f123","#f5a527","#db2f14")) +
scale_fill_manual(values=c("#3aba1e","#a1da27","#f0f123","#f5a527","#db2f14")) +
theme_bw() +
theme(strip.background=element_blank()) +
labs(x ="Grazing pressure", y ="Species diversity (richness)") +
facet_grid(~Location)
fig_eshannon <- ggplot(aes(x = Plot_rev, y = EShannon, fill = Plot_rev), data = Plot_variables_s) +
#geom_violin(trim = FALSE, alpha = .25, colour = NA) +
geom_jitter(width = .15, aes(color = Plot_rev), alpha = .5, shape = 16, size = 3) +
stat_summary(fun.data = mean_cl_boot, geom = "errorbar", width = .15, color = "grey25") +
stat_summary(fun.y = mean, geom = "line", size = .5, aes(group = Location), color = "grey25") +
stat_summary(fun = mean, geom = "point", size = 4, aes(fill = Plot_rev), shape = 22) +
scale_color_manual(values=c("#3aba1e","#a1da27","#f0f123","#f5a527","#db2f14")) +
scale_fill_manual(values=c("#3aba1e","#a1da27","#f0f123","#f5a527","#db2f14")) +
theme_bw() +
theme(strip.background=element_blank()) +
labs(x ="Grazing pressure", y ="Species diversity (hill numbers)") +
facet_grid(~Location)
#Print
grid.arrange(fig_rich, fig_eshannon, nrow = 2)
fig_pH <- ggplot(aes(x = Plot_rev, y = pH, fill = Plot_rev), data = Plot_variables_s) +
#geom_violin(trim = FALSE, alpha = .25, colour = NA) +
geom_jitter(width = .15, aes(color = Plot_rev), alpha = .5, shape = 16, size = 3) +
stat_summary(fun.data = mean_cl_boot, geom = "errorbar", width = .15, color = "grey25") +
stat_summary(fun.y = mean, geom = "line", size = .5, aes(group = Location), color = "grey25") +
stat_summary(fun = mean, geom = "point", size = 4, aes(fill = Plot_rev), shape = 22) +
scale_color_manual(values=c("#3aba1e","#a1da27","#f0f123","#f5a527","#db2f14")) +
scale_fill_manual(values=c("#3aba1e","#a1da27","#f0f123","#f5a527","#db2f14")) +
theme_bw() +
theme(strip.background=element_blank()) +
labs(x ="Grazing pressure", y ="Soil fertility (pH)") +
facet_grid(~Location)
fig_EC <- ggplot(aes(x = Plot_rev, y = EC, fill = Plot_rev), data = Plot_variables_s) +
#geom_violin(trim = FALSE, alpha = .25, colour = NA) +
geom_jitter(width = .15, aes(color = Plot_rev), alpha = .5, shape = 16, size = 3) +
stat_summary(fun.data = mean_cl_boot, geom = "errorbar", width = .15, color = "grey25") +
stat_summary(fun.y = mean, geom = "line", size = .5, aes(group = Location), color = "grey25") +
stat_summary(fun = mean, geom = "point", size = 4, aes(fill = Plot_rev), shape = 22) +
scale_color_manual(values=c("#3aba1e","#a1da27","#f0f123","#f5a527","#db2f14")) +
scale_fill_manual(values=c("#3aba1e","#a1da27","#f0f123","#f5a527","#db2f14")) +
theme_bw() +
theme(strip.background=element_blank()) +
labs(x ="Grazing pressure", y ="Soil fertility (EC)") +
facet_grid(~Location)
fig_OM <- ggplot(aes(x = Plot_rev, y = MO, fill = Plot_rev), data = Plot_variables_s) +
#geom_violin(trim = FALSE, alpha = .25, colour = NA) +
geom_jitter(width = .15, aes(color = Plot_rev), alpha = .5, shape = 16, size = 3) +
stat_summary(fun.data = mean_cl_boot, geom = "errorbar", width = .15, color = "grey25") +
stat_summary(fun.y = mean, geom = "line", size = .5, aes(group = Location), color = "grey25") +
stat_summary(fun = mean, geom = "point", size = 4, aes(fill = Plot_rev), shape = 22) +
scale_color_manual(values=c("#3aba1e","#a1da27","#f0f123","#f5a527","#db2f14")) +
scale_fill_manual(values=c("#3aba1e","#a1da27","#f0f123","#f5a527","#db2f14")) +
theme_bw() +
theme(strip.background=element_blank()) +
labs(x ="Grazing pressure", y ="Soil fertility (OM)") +
facet_grid(~Location)
fig_C <- ggplot(aes(x = Plot_rev, y = TC, fill = Plot_rev), data = Plot_variables_s) +
#geom_violin(trim = FALSE, alpha = .25, colour = NA) +
geom_jitter(width = .15, aes(color = Plot_rev), alpha = .5, shape = 16, size = 3) +
stat_summary(fun.data = mean_cl_boot, geom = "errorbar", width = .15, color = "grey25") +
stat_summary(fun.y = mean, geom = "line", size = .5, aes(group = Location), color = "grey25") +
stat_summary(fun = mean, geom = "point", size = 4, aes(fill = Plot_rev), shape = 22) +
scale_color_manual(values=c("#3aba1e","#a1da27","#f0f123","#f5a527","#db2f14")) +
scale_fill_manual(values=c("#3aba1e","#a1da27","#f0f123","#f5a527","#db2f14")) +
theme_bw() +
theme(strip.background=element_blank()) +
labs(x ="Grazing pressure", y ="Soil fertility (TC)") +
facet_grid(~Location)
fig_N <- ggplot(aes(x = Plot_rev, y = TN, fill = Plot_rev), data = Plot_variables_s) +
#geom_violin(trim = FALSE, alpha = .25, colour = NA) +
geom_jitter(width = .15, aes(color = Plot_rev), alpha = .5, shape = 16, size = 3) +
stat_summary(fun.data = mean_cl_boot, geom = "errorbar", width = .15, color = "grey25") +
stat_summary(fun.y = mean, geom = "line", size = .5, aes(group = Location), color = "grey25") +
stat_summary(fun = mean, geom = "point", size = 4, aes(fill = Plot_rev), shape = 22) +
scale_color_manual(values=c("#3aba1e","#a1da27","#f0f123","#f5a527","#db2f14")) +
scale_fill_manual(values=c("#3aba1e","#a1da27","#f0f123","#f5a527","#db2f14")) +
theme_bw() +
theme(strip.background=element_blank()) +
labs(x ="Grazing pressure", y ="Soil fertility (TN)") +
facet_grid(~Location)
fig_CN <- ggplot(aes(x = Plot_rev, y = CN, fill = Plot_rev), data = Plot_variables_s) +
#geom_violin(trim = FALSE, alpha = .25, colour = NA) +
geom_jitter(width = .15, aes(color = Plot_rev), alpha = .5, shape = 16, size = 3) +
stat_summary(fun.data = mean_cl_boot, geom = "errorbar", width = .15, color = "grey25") +
stat_summary(fun.y = mean, geom = "line", size = .5, aes(group = Location), color = "grey25") +
stat_summary(fun = mean, geom = "point", size = 4, aes(fill = Plot_rev), shape = 22) +
scale_color_manual(values=c("#3aba1e","#a1da27","#f0f123","#f5a527","#db2f14")) +
scale_fill_manual(values=c("#3aba1e","#a1da27","#f0f123","#f5a527","#db2f14")) +
ylim(0,40) +
theme_bw() +
theme(strip.background=element_blank()) +
labs(x ="Grazing pressure", y ="Soil fertility (CN)") +
facet_grid(~Location)
fig_P <- ggplot(aes(x = Plot_rev, y = P_Olsen, fill = Plot_rev), data = Plot_variables_s) +
#geom_violin(trim = FALSE, alpha = .25, colour = NA) +
geom_jitter(width = .15, aes(color = Plot_rev), alpha = .5, shape = 16, size = 3) +
stat_summary(fun.data = mean_cl_boot, geom = "errorbar", width = .15, color = "grey25") +
stat_summary(fun.y = mean, geom = "line", size = .5, aes(group = Location), color = "grey25") +
stat_summary(fun = mean, geom = "point", size = 4, aes(fill = Plot_rev), shape = 22) +
scale_color_manual(values=c("#3aba1e","#a1da27","#f0f123","#f5a527","#db2f14")) +
scale_fill_manual(values=c("#3aba1e","#a1da27","#f0f123","#f5a527","#db2f14")) +
theme_bw() +
theme(strip.background=element_blank()) +
labs(x ="Grazing pressure", y ="Soil fertility (Available P)") +
facet_grid(~Location)
fig_NP <- ggplot(aes(x = Plot_rev, y = NP, fill = Plot_rev), data = Plot_variables_s) +
#geom_violin(trim = FALSE, alpha = .25, colour = NA) +
geom_jitter(width = .15, aes(color = Plot_rev), alpha = .5, shape = 16, size = 3) +
stat_summary(fun.data = mean_cl_boot, geom = "errorbar", width = .15, color = "grey25") +
stat_summary(fun.y = mean, geom = "line", size = .5, aes(group = Location), color = "grey25") +
stat_summary(fun = mean, geom = "point", size = 4, aes(fill = Plot_rev), shape = 22) +
scale_color_manual(values=c("#3aba1e","#a1da27","#f0f123","#f5a527","#db2f14")) +
scale_fill_manual(values=c("#3aba1e","#a1da27","#f0f123","#f5a527","#db2f14")) +
theme_bw() +
theme(strip.background=element_blank()) +
labs(x ="Grazing pressure", y ="Soil fertility (ln(NP))") +
facet_grid(~Location)
#Print
grid.arrange(fig_pH, fig_EC, fig_OM, fig_C, fig_N, fig_CN, fig_P, fig_NP, nrow = 8)